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The basic equations are derived for compressible flow in a stepped labyrinth gas seal. 
The flow is assumed to be completely turbulent in the circumferential direction where 
the friction factor is determined by the Blasius relation. Linearized zeroth and first-order 
perturbation equations are developed for small motion about a centered position by an 
expansion in the eccentricity ratio. The zeroth-order pressure distribution is found by 
satisfying the leakage equation while the circumferential velocity distribution is determined 
by satisfying the momentum equations. The first order equations are solved by a separation 
of variables solution. Integration of the resultant pressure distribution along and around 
the seal defines the reaction force developed by the seal and the corresponding dynamic 
coefficients. 

The results of this analysis are presented in the form of a parametric study, since there 
are no known experimental data for the rotordynamic coefficients of stepped labyrinth gas 
seals. The parametric study investigates the relative rotordynamic stability of convergent, 
straight and divergent stepped labyrinth gas seals. The results show that, generally, the 
divergent seal is more stable, rotordynamically, than the straight or convergent seals. The 
results also show that the teeth-on-stator seals are not always more stable, rotordynam- 
ically, than the teeth-on-rotor seals as was shown by experiment by Childs and Scharrer 
(1986b) for a 15 tooth seal. 


INTRODUCTION 

The problem of self excited vibration in turbomachinery due to labyrinth seals has led to the 
development of many analyses for the straight-through type of labyrinth seal, e.g. Childs 
and Scharrer (1986a), Iwatsubo (1980), Jenny et al. (1984) and Scharrer (1987). There 
has also been an extensive amount of experimental data collected for the rotordynamic 
coefficients of straight-through labyrinth seals, e.g. Childs and Scharrer (1986b), Leong 
(1983), Scharrer (1987) and Wachter and Benckert (1980). The stepped labyrinth seal is 
as prevalent in turbomachinery as the straight-through labyrinth seal; however, there has 
been very little or no experimental or analytical work done to quantify the rotordynamic 
characteristics of this type of seal. This paper presents an analysis of the stepped labyrinth 
seal which is a derivative of the straight-through analysis of Childs and Scharrer (1986a). 
The results of this analysis will be presented in the form of a parametric study, since there 
are no known experimental data. 
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PROCEDURE 


The analysis presented here is developed for the stepped type of labyrinth seal shown 
in figure 1. The continuity and momentum equations will be derived for a single cavity 
control volume as shown in figures 2,3,4 and 5. A leakage model will be employed to 
account for the axial leakage. The governing equations will be linearized using perturbation 
analysis for small motion about a centered position. The zeroth-order continuity and 
momentum equations will be solved to determine the steady state pressure and velocity 
for each cavity. The first-order continuity and momentum equations will be reduced to 
linearly independent, algebraic equations by assuming an elliptical orbit for the shaft 
and a corresponding harmonic response for the pressure and velocity perturbations. The 
force and force coefficients for the seal are found by integration of the first-order pressure 
perturbation along and around the shaft. 


ASSUMPTIONS 


(1) The fluid is an ideal gas. 

(2) Pressure variations within a chamber are small compared to the pressure differences 
across a seal strip. 

(3) The lowest frequency of acoustic resonance in the cavity is much higher than that 
of the rotor speed. 

(4) The eccentricity of the rotor is small compared to the radial seal clearance. 

(5) Although the shear stress is significant in the determination of the flow parameters 
(velocity, etc.), the shear stress forces on the rotor are small when compared to the pressure 
forces. 

(6) The cavity flow is turbulent and isothermal. 

(7) Added mass terms are neglected. 

(8) The seal strip is midway between steps. 

GOVERNING EQUATIONS 


Continuity Equation 

The control volumes of figures 2 and 3 have a unit circumferential width. Their continuity 
equation is 


8{pA) 1 8(pVA) 

dt Rs{ dO 


Rfi + 1 
+ Rsi 


m^ + i — th{ = 0 


where the transverse surface area, A, is defined by 


Ai — ( Bi + Hi + Bi + i + Hi+i)L{/2 


( 1 ) 


Momentum Equations 

The momentum equation (2) is derived using figures 4 and 5 which show the pressure 
forces and shear stresses acting on the control volume. 


178 



dpViAi 2pViAi dVi pV^dA ViAdp Rsj^thj^ 

dt + Rsi dO + Rsi dO + Rsi d0 + Rsi 

A. Qp. 

— ihiVi - 1 = + T riO-U{Li + <0 — T t iaSi(Li + d) (2) 

The shear stress model is the same as that used by Childs and Scharrer (1986a). The 
differences are in the dimensionless shear stress lengths and the hydraulic diameter. The 
dimensionless shear stress length is defined for the teeth-on-rotor labyrinth by 


oSj = 1 ari = (2J3; + L i + d)/(Li + <f) 

(3) 

and for teeth-on-stator seals by 


ari = 1 asi = (2 Bi + + d)/(Li + d) 

(4) 

The hydraulic diameter, is defined by 


2 (Bi + Hi){Li + d) 
h% Hi+Bi + Li + d 

(5) 

Reduced Equations 

If equation (1) times the circumferential velocity is now subtracted from equation (2), the 
following reduced form of the momentum equation is obtained: 


dVi pViAi dVi 

+ 


dt ' Rsi 
Ai dPi 


Rsi dO 


dO 

+ TridTiLi 


+ +m,(Vi-V-_ 1 ) 

TgidSiLi 


( 6 ) 


The number of variables is reduced by using the ideal gas law to eliminate the density 
terms. 


Pi = PiRT 


( 7 ) 


Leakage Equation 

The leakage equation (8), flow coefficient, kinetic energy carryover coefficient, and the 
choked flow solution algorithm are the same as those used by Childs and Scharrer (1986a). 
The kinetic energy carryover coefficient has the added restriction of being unity when the 
step height is greater than the clearance. 


rhi = puPiiHi 



( 8 ) 


Perturbation Analysis 

The solution procedure will be summarized here as the details are given in Childs and 
Scharrer (1986a). Introduction of the following perturbation variables into the governing 
equations yields zeroth and first order continuity and momentum equations. 
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Pi = Pai + ePu Hi = Cu + eHi 

Vi = Voi + eV u Ai = A 0 + eLiH t 

where e = e 0 /Cri is the eccentricity ratio. 

Zeroth- Order Solution 

The zeroth-order leakage equation is 


Rs i+l rhi + 1 _ ^ ^ 

— — TYli — T7l 0 

RS{ 

and can be solved iteratively to yield zeroth-order cavity pressure values. 
The zeroth-order circumferential momentum equation is 


( 9 ) 


*o{RH+lV ei /RSi loi-l) — ( TroiQTi T to i<lSi)Li (10) 

and is solved using a Newton root finding technique to yield zeroth-order cavity circum- 
ferential velocity components. 


First- Order Solution 

The governing first-order equations (11,12), define the pressure and velocity fluctuations 
resulting from the seal clearance function. The continuity and momentum equations follow 
in order: 


G u 


dP 


li 


at 


+ Gu 


Voi dP U 

Rsi dO 


+ 


R7,~ar + GziPu 


+ + GsiPii+1 = —GoiHu — <?2i — <?2i‘ °* U 


at 


Rsi ao 


v dVu XuVoi dVu Aoi dPu Tr ^ 

Xu ~m R^~df + RT i ~dF + x * Vu ~ 

+ XsiPu + XuPli-l = XsiHu 


(11) 


( 12 ) 


where the X’s and G’s are defined in Appendix A. With an assumption of an elliptical shaft 
orbit, these equations can be reduced to the following system of linear algebraic equations: 

IA-ilpfi -0 + ww + lA+i](x i+1 ) = | (ft) + 1 (CO ( 13 ) 

where 


(Xi) = (P +, , Pi , p - , p - , v+ , v + , v^,7 , V - ) T 


)T 

) T 


The A matrices and column vectors B and C are given in Childs and Scharrer (1986a). To 
use equation (13) for the entire solution, a system matrix can be formed which is block 
tridiagonal in the A matrices. The size of this resultant matrix is (8NC X 8NC) since 
pressure and velocity perturbations at the inlet and the exit are assumed to be zero. This 
system is easily solved by various linear equation algorithms, and yields a solution of the 
form: 
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p+ = -F + - + -Ff- 

* st * ast ' ^ oat 

a _ 6 __ 

-P« = + ~ F bsi 

p+ = -F + . + - F& 

* ct ^ * act 1 ^ * bet 

a 6 

■Pet = ~ F ad + “-Pfect 


(14) 


DETERMINATION OF DYNAMIC COEFFICIENT 

The force-motion equations for a labyrinth seal are assumed to be of the form 

-{£}-[-" ;]{*} 

The solution of equation (15) for the stiffness and damping coefficients is the objective 
of the current analysis. The solution procedure used for this analysis is the same one 
used by Childs and Scharrer (1986a). The desired solution for the stiffness and damping 
coefficients is 


NG 


t=l 

NG 

t=l 

E> NC 

C = -—J2(F^ i + F~ i )L i 


c = 


w 

7tRs 

u 


t=l 

NG 


£(*& + w 


t=l 


(16) 


RESULTS 

The geometry used in the parametric study is given in table 1. The pitch of the teeth, 
the step height, the radial clearance, the tooth height and the inlet radius were all kept 
constant. The number of teeth were varied between 5 and 15 for the three types of 
seals investigated: straight, converging and diverging which are shown in figure 6. The 
operating conditions used for the study are given in table 2. The only variable is the inlet 
circumferential velocity ratio which was varied between 0.25 and 1.0. 
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Table 1. Basic Geometry Studied 


Radial Clearance, Cr 
Tooth Height, B 
Tooth Pitch, L 
Seal Radius, Rs 
Tooth Width, tp 
Step H eight, d 
Number of 


= 0.127 mm (0.005 in) 

= 3.175 mm (0.125 in) 

= 2.175 mm (0.125 in) 

= 75.6 mm (2.979 in) 
= 0.152 mm (0.006 w) 

= 1.000 mm (0.04 in) 

Teeth, NT = 5, 10, 15 


Table 2. Operating Conditions 

Reservoir Pressure, Pr = 7.0 bar (101.0 psi) 
Sump Pressure, Ps = 1.01 bar (14.7 psi) 
Shaft Speed,u> = 20,000 cpm (rpm) 
Temperature, T = 300 K (80 deg F) 

Inlet Swirl Ratio,Vi/Rsw = 0.25, 0.5, 0.75, 1.0 

Fluid is air 


Leakage 

Figures 7 and 8 show leakage versus number of teeth and configuration for teeth-on-rotor 
and teeth-on-stator seals, respectively. The figures show that the leakage decreases with 
increasing number of teeth for both teeth-on-rotor and teeth-on-stator seals. The figures 
also show that the diverging seal yields maximum leakage while the converging seal has 
the minimum leakage. This result is expected since the effective area of the diverging seal 
is greater than that of the converging or straight seal. 

Stiffness 

Figures 9 and 10 show direct stiffness versus number of teeth and configuration for an inlet 
circumferential velocity ratio of 1.0 for teeth-on-rotor and teeth-on-stator seals, respec- 
tively. Figure 9 shows that for converging and straight teeth-on-rotor seals a maximum 
stiffness occurs for the 10 tooth seal. For the teeth-on-stator seal, figure 10 shows that 
there is an optimum only for the diverging seal. A maximum may occur for geometries 
other than those investigated in this study. This phenomenon did not occur for inlet cir- 
cumferential velocity ratios of 0.25 or 0.5. Both figures also show that the diverging seal 
yielded the largest value of stiffness. This was true for all of the geometries and operating 
conditions investigated. 

Figures 11 and 12 show cross-coupled stiffness versus configuration and number of teeth for 
an inlet circumferential velocity ratio of 1.0 for teeth-on-rotor and teeth-on-stator seals, 
respectively. Figure 11 shows that the straight teeth-on-rotor seal yields the maximum 
cross-coupled stiffness for both the 10 and 15 tooth seals. This phenomenon did not occur 
for the 0.25 inlet circumferential velocity ratio case. Figure 12 shows that the same is 
true for the 10 tooth teeth-on-stator seal. This phenomenon did not occur for any inlet 
circumferential velocity ratio case other than 1.0. The figures also show that the cross- 
coupled stiffness increases as the number of teeth increase. 
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Damping 

Figures 13 and 14 show direct damping versus configuration and number of teeth for 
an inlet circumferential velocity ratio of 1.0 for teeth-on-rotor and teeth-on-stator seals, 
respectively. Figure 13 shows that the straight teeth-on-rotor seal yields a maximum value 
of damping for the 10 and 15 tooth seals. This did not occur for an inlet circumferential 
velocity ratio of 0.25 or 0.5. Figure 14 shows that direct damping is a maximum for the 
diverging teeth-on-stator seal and a minimum for the converging seal for the 10 and 15 
tooth seals. The opposite is true for the 5 tooth seal. The figures also show that direct 
damping increases as the number of teeth increases. 

Relative Rotordynamic Stability 

The relative stability of the seal configurations investigated will be compared using the 
whirl frequency ratio which is a ratio of the destabilizing influence divided by the stabilizing 
influence and is defined as 


Jl 

<jJC 

When this ratio is greater than 1.0, the seal is a destabilizing influence. The seal with 
the lowest value is considered to be the most stable. Figures 15 and 16 show the whirl 
frequency ratio versus seal configuration for inlet circumferential velocity ratios of 0.25 
and 1.0, respectively. Both figures show that the teeth-on-rotor seal is more stable than 
the teeth-on-stator seal for the 5 tooth configuration while the opposite is true for the 15 
tooth configuration. The latter result is supported by the test data of Childs and Scharrer 
(1986b). The 10 tooth configuration seems to be in a transition region. The figures also 
show that except for the 5 tooth seals at inlet swirls of 1.0, the diverging seal configuration 
is more stable than the converging and straight configurations. A more complete stability 
analysis would show that the converging seal may in some cases be a more unstable seal 
due to the large negative direct stiffness values illustrated in figures 9 and 10. 


CONCLUSIONS 

This paper has presented an analysis for the rotordynamic coefficients of stepped labyrinth 
gas seals. The results of this analysis were presented in the form of a parametric study, 
since no known experimental data existed. The results of this study support the following 
conclusions: 

1) The converging stepped labyrinth seal leaks less than either the straight-through 
or the divergent seals. 

2) The diverging stepped labyrinth seal had a higher value of direct stiffness than either 
the straight-through or convergent seals for both the teeth-on-rotor and teeth-on-stator 
case. 

3) A maximum direct stiffness was obtained for a 10 tooth convergent teeth-on-rotor 
seal and a 10 tooth divergent teeth-on-stator seal for inlet circumferential velocity ratios 
of 0.75 or larger. 

4) Cross-coupled stiffness increases as the number of teeth, and therefore length (in 
this study), increase. 
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5) Direct damping increases as the number of teeth, and therefore length (in this 
study), increase. 

6) For most of the cases studied, the diverging seal is more stable, rotordynamically, 
than either the straight-through or converging seals for both teeth-on-rotor and teeth-on- 
stator seals. 

7) For the 5 tooth seals, the teeth-on-rotor seal was more stable, rotordynamically, 
than the teeth-on-stator seals. The opposite was true for the 15 tooth seals. 


NOMENCLATURE 


A Cross-sectional area of control volume (L 2 ); illustrated in figure (3) 

B Height of labyrinth seal strip (L); illustrated in figure (1) 

C Direct damping coefficient (Ft/L) 

Cr Nominal radial clearance (L); illustrated in figure (1) 

Dh Hydraulic diameter of cavity (L); introduced in equation (5) 

H Local radial clearance (L) 

K Direct stiffness coefficient (F/L) 

L Pitch of seal strips (L); illustrated in figure (1) 

NT Number of seal strips 

NC = NT — 1 Number of cavities 

P Pressure [F/L 2 ) 

R Gas constant [L 2 /Tt 2 ) 

Rs Radius of control volume (L); illustrated in figure (l) 

Rsu Surface velocity of rotor (L/t) 

T Temperature (T) 

Tp Tooth tip width (L); illustrated in figure (l) 

V Average circumferential velocity for control volume (L/t); 

illustrated in figure (2) 

or, as Dimensionless length upon which shear stress acts; introduced 

in equation (3) and (4) 

c Cross coupled damping coefficient (Ft/L); in equation (16) 

d Step height (L) , illustrated in figure 1 

e a Displacement of the seal rotor from centered position (L) 

k Cross coupled stiffness coefficient (F/L); in equation (16) 

m Leakage mass flow rate per circumferential length (M/Lt) 

mr, nr, ms, ns Coefficients for friction factor 

t Time (t) 

(jj Shaft angular velocity (l/t) 

p Density of fluid ( M/L 3 ) 

v Kinematic viscosity [L 2 /t) 

€ = e 0 /Cr Eccentricity ratio 
7 Ratio of specific heats 


C-3 
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Subscripts 


0 Zeroth-order component 

1 First-order component 

i i-th chamber value 

x X-direction 

y Y-direction 

r Reservoir value 

s Sump value 


APPENDIX A: DEFINITION OF FIRST-ORDER COEFFICIENTS 


G, = — ; G 2 = 

1 RT 2 RT 


G* = 


rh 0 P 0 


3 p2 __ p2 

r oi - 1 r oi 

rn 0 pu{l -1) 


rh 0 pu+i{l- l), c 4C , \( p oi\ 1 1 

„ ^p~ + . (5 - 4Sii+i) l^rj “sr 


+ 


, ''MCin I •‘//p 4C , \/c< WoPoi , ihoPoi Rsi+l 

+ ™ T : (5 “ + 1 )+ p2 p2 + p2 _ p2 p— 

“ni r oi+l ^ai r oi + 1 


G 4 = 


G5 = 


-rhoPvi - 1 rh 0 pu (7 - 




oi * oi+1 * 01 * oi+l 

-=JL 
1 


Pi- 1 - -P5 ’H p « 

-m«,P CT -n Xsi+i - 1) 


^-■P3 +1 




G 6 =m, 


Gr» - Cry_i 

’ CnCn-i 


x, = 


i2T 


x, = - ' h ° ( pt ^I-Oh 1 ) p li (4S li - 5)(S U + 1) 

“oi— 1 'oi K~l“oi 


X 4 = 


r T i ar»Lj(l + mr) i r„(l + ms)Li as» 

^ + Poi 

m 0 (^ - — Vpt- 1 )(7 ~ l) 






FU^Su-S) (^i) 


v _ . r^asjLi(2 + ms) r r jarjLt(2 4- mr) 

^2 m© "l I - / JTJ T/ \ 

Yoi (Ps»w - V*) 

y- w*o(loi Fo»— i ) TeiasiLimsDhi T r j(iriLitnT D m 

5 ~ 2(Gr< + Bi) 2 2(Cr t + J?<) 2 
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APPENDIX B: SEPARATION OF GOVERNING EQUATIONS 

Continuity: 

cos(0 + ut) :G l PZ (u. + VJ + G,Pt + G.P+., 

+ G 5 P+ +X = Go^-^ 

sin(0 + wt) : — GiP^ fw + ^ — G x + GsP.* + G 4 P„;_ 1 

+ G ^ = f® + 

cos(0 — wt) : - GiP,7 ^ G ^-i 

(0 + 6 ) 

+ G 5 P c J +l =G 6 i-j-t 

sin(0 -ut):- GtP- (« - j£) - G,^ KJ + G,P,7 + G 4 J5_, 

+ GsP.7 +1 = f (-£)(• + ») 


Momentum: 


cos(0 + wt) :XiV^ + w^ + -j^-P,"i + -^^cT + ^3^c* 


+ X 4 P+_ 1 = X 5 


(6 -a) 


m(0 + wt) : — XiV£ rh 0 V ^_ 1 + X3 P ti 

+ X 4 P.t-i = 0 

os(0 — vt) :XiV,7 _ + l^s-P ** XzVd ~ WoVci - 1 + XqPqi 


+ X 4 P C ^_ 1 — X 5 


(a + 6) 


«n(< - art) :X:Vu (« - W - -M'.t., + Xtfi 


+ X 4 P 8 i_! — 0 
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DEFINITION OF MATRIX ELEMENTS 
Ai-t MATRIX 

-A 1,2 = A 2 , l = As, 4 = A4,3 = G 4 

A5,2 = A 6 ,l = A 7,4 = As, s = X 4 
A5,6 = A 6 ,5 = A 7 ,s = As, 7 = -n»o 
The remaining elements are zero. 

Ai MATRIX 



Al ,2 = A 2 ,i = A 3 , 4 = A 4 ,3 = Gs 
As, 2 = A 6 ,i = A 7 ,4 = As, 3 = X3 

A54 = -As, 2 = A 7 ,3 = -As, 4 = ^7 



A5 , 6 = As, 5 = A 7 ,8 — As, 7 — X2 


The remaining elements are zero. 

Ai + 1 MATRIX 
Al,2 = A24 = A3, 4 = A4,3 = Gs 
The remaining elements are zero. 

B and C Column Vectors 
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Figure 1. A Typical Stepped Labyrinth Seal. 
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Figure 4. Forces on Control Volume. 




Figure 5. Pressure Forces on Control Volume. 
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Figure 6. Configurations Studied. 
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STEPPED LABYRINTH SEAL 
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Figure 7. Leakage versus number of teeth and seal configuration for a 
teeth-on-rotor seal. 


STEPPED LABYRINTH SEAL 


TEETH— ON- ST AT OR 



NUMtlCR OP TEETH _ 

O CONVERCINO ♦ SIRAICHT • OIVERCINC 

Figure 8. Leakage versus number of teeth and seal configuration for a 
teeth-on-stator seal. 
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STEPPED LABYRINTH SEAL 



NUMBER or TEETH 

O CONVERClNC ■» STRAICHT ♦ DfVERCINC 

Figure 9. Direct Stiffness versus number of teeth and seal configuration 
for a teeth-on- rotor seal. Inlet circumferential velocity ratio is 1.0. 


STEPPED LABYRINTH SEAL 



NUMEER cr TEETH 

D CONVERClNC * STRAICHT « DFVERCINC 

Figure 10. Direct Stiffness versus number of teeth and seal configuration 
for a teeth-on-stator seal. Inlet circumferential velocity ratio is 1.0. 
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STEPPED LABYRINTH SEAL 



CONFIGURATION 

D 3 TEETH ♦ 10 TEETH • 13 TEETH 

Figure 11. Cross-coupled stiffness versus seal configuration and number 
of teeth for a teeth-on-rotor seal. Inlet circumferential velocity ratio is 
1 . 0 . 


STEPPED LABYRINTH SEAL 

TEETH-ON-STATOR V/RS- - 1.0 



CONFIGURATION 
10 TCETH 


3 TEETH 4 10 TCETH • 13 TEETH 

Figure 12. Cross-coupled stiffness versus seal configuration and number 
of teeth for a teeth-on-stator seal. Inlet circumferential velocity ratio is 
1 . 0 . 
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Figure 13. Direct damping versus seal configuration and number of teeth 
for a teeth-on-rotor seal. Inlet circumferential velocity ratio is 1.0. 


STEPPED LABYRINTH SEAL 



. _ CONFIGURATION 

O 5 TEETH ♦ 10 TEETH # 13 TEETH 

Figure 14. Direct damping versus seal configuration and number of teeth 
for a teeth-on-rotor seal. Inlet circumferential velocity ratio is 1.0. 
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STEPPED LABYRINTH SEAL 



Figure 15. Whirl frequency ratio versus seal configuration for an inlet 
circumferential velocity ratio of 0.25. 


STEPPED LABYRINTH SEAL 
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Figure 16. Whirl frequency ratio versus seal configuration for an inlet 
circumferential velocity ratio of 1.0. 
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